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1. INTRODUCTION 

The intensity and amplitude statistics of SAR images, such as L-Band HH for 
SEASAT and S1R-B, and C-Band VV for ERS-1 have been extensively investigated 
for various terrain, ground cover and ocean surfaces. Less well-known are the 
statistics between multiple channels of polarimetric or interferometric SARs, 
especially for the multi-look processed data. In this paper, we investigate the 
probability density functions (PDFs) of phase differences, the magnitude of complex 
products and the amplitude ratios, between polarization channels (i.e. HH,HV, and 
VV) using 1-look and 4-look AIRSAR polarimetric data. Measured histograms are 
compared with theoretical PDFs which were recently derived based on a complex 
Gaussian model (Lee et al., 1993). 

NASA/JPL 1-look and 4-look AIRSAR data of Howland Forest and San 
Francisco were used for comparison. Histograms from 1-look SAR data agreed with 
theoretical PDFs. However, discrepancies were found when matching the 4-look 
polarimetric data with the 4-look PDFs. Instead, We found that the 3-look PDFs 
matched better. The problem was traced to the averaging of correlated 1-look pixels. 
We also verified these theoretical PDFs for forest, ocean surfaces, park areas and city 
blocks. This indicates that ground cover and terrain with different scattering 
mechanisms can be represented with this statistical model. 

2. THE COMPLEX GAUSSIAN MODEL 

A polarimetric radar measures the complete scattering matrix S of a medium 
at a given incidence angle. For a reciprocal medium, the three unique complex 
elements are S hh , S hv , and S^. Circular Gaussian Conditions (Goodman, 1985) are 
assumed. The circular Gaussian assumption has been verified by Sarabandi(I992) 
using 1-look polarimetric SAR data. 

The multi-look AIRSAR processor compresses polarimetric data by averaging 
Mueller matrices of 1 -look pixels in the azimuth direction. Due to oversampling, the 
neighboring 1-look pixels of AIRSAR data are somewhat correlated. For mathematical 
simplicity, the PDFs for multi-look phase difference, etc. were derived under the 
assumption of statistical independence. 

3. MULTI-LOOK PHASE DIFFERENCE DISTRIBUTION 

Let n be the number of looks, and S^k) and S 2 (k) be the kth 1-look samples 
of any two components of the scattering matrix. For polarimetric and interferometric 
radars, the multi-look phase difference is obtained by 


33 



( 1 ) 


y=Aig[± £ s x (k) s' 2 {k) ] 

n k=l 


Under the circular Gaussian assumption, the multi-look phase distribution 
were derived by multiple integrations of special functions (Lee et al., 1993). The 
multi-look phase difference PDF is 
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where F(n, 1; 1/2; 0 2 ) is a Gauss hypergeometric function, and 6 and |p c | are the phase 
and the magnitude of the complex correlation coefficient defined as 
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The 1-look (n=l) PDF can be obtained by applying mathematical identities of 
the hypergeometrical function (Lee et al., 1 993). The 1-look PDF obtained is identical 
to that of Kong (1988) and Sarabandi ( 1992). The PDF of Eq. (2) depends only on the 
number of looks and the complex correlation coefficient. The peak of the distribution 
is located at 0=0. A plot of standard deviation versus |p c | is given in Fig. 1, which 
verifies a similar but less accurate figure (Zebker, 1992) computed with 
interferometric radar data. As shown in Fig. I, multi-look processing effectively 
reduces the phase error, especially when n=16 and 32. 

4. DISTRIBUTION OFTHE MAGNITUDE OFTHE NORMALIZED MULTI-LOOK 
COMPLEX PRODUCT 

The magnitude of product of S x and S 2 is an important measure in 
polarimetric SAR, and it is the magnitude of interferogram from an interferometric 
SAR. The normalized magnitude is defined as 

n 

l~E S^k) S* 2 {k) | 

' (5) 

v/stlsj 2 ] £t|s 2 | 2 ] 


The PDF of £ (Lee et al., 1 993) is 
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where / 0 ( ) and K n ( ) are modified Bessel functions. The PDF for the unnormalized 
magnitude can be easily obtained from Eq.(6) by using Eq. (5). The standard 


34 



deviations of this PDF is plotted versus the correlation coefficient in Fig. 2. This plot 
indicates that the magnitude of complex product can be used as a discriminator, 
especially for large n. 

5. MULTI-LOOK AMPLITUDE RATIO DISTRIBUTIONS 

The amplitude ratio between S hh and S vv has been an important discriminator 
in the study of polarimetric radar returns. Let the normalized ratio 
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The multi-look amplitude ratio PDF (Lee et al., 1993) is 
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For n=l, we have the 1-look amplitude distribution, which is identical to the results 
of Kong(1988), when his PDF is properly normalized. 



Fig. 1 Phase difference standard deviations 
versus the correlation coefficient, tp c U for 
the number of look from n= 1, 2, 4, 8, 16 
and 32. The value of 0 = 0 . 
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Fig. 2 Standard deviations of Normalized 
magnitude of multi-look interferogram 
versus correlation coefficient, |p c |, for the 
number of look from n= 1 , 2, 4, 8 and 16. 


6. COMPARISON OF THEORETICAL PDFs WITH AIRSAR MEASUREMENTS 

In this section, the AIRSAR polarimetric data was used to compute histograms 
of phase differences, normalized products and amplitude ratios to verify this complex 
Gaussian model. Homogeneous regions of forest were selected, and the complex 
correlation coefficient (i.e., 6 and |p c |) and histograms were computed. We have 
checked the 1 -look C-Band and L-Band data, and found that the agreements between 
histograms and their corresponding PDFs are good. However, discrepancies exist in 
the 4-look data. Using 4-look Howland Forest data (CM1084), the match between the 
histogram and the 4-look phase difference PDF are not as good (Fig. 3A). A better 
match was found with a 3-look PDF (Fig. 3B). Histograms of amplitude and 
normalized product are also shown good agreements with 3-look PDFs. To save space, 
only the case of HH to VV ratio is shown in Fig. 3C. 
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The problem was traced to the averaging of correlated 1-look neighboring 
pixels during the multi-look processing. To verify it, we use 1-look data of Howland 
Forest (HR1084C) and perform the 4-look processing by averaging four pixels 
separated by two pixels in the azimuth direction. Since the correlation between every 
other pixels is much less than that between its immediate neighbors, statistical 
independence is assured. The results are shown in Fig.4 for all three variables under 
study. The agreement with 4-look PDFs is very good. We can conclude that due to 
the correlation of 1-look data, the 4-look A1RSAR data has the characteristics of a 
3 -look. 



Fig. 3 Experimental results using 4-look A1RSAR data of Howland Forest (CM1804) 
with |p c |=0.491 and 0=84.9°. (A) The histogram of phase difference between HH and 
VV and the theoretical 4-look PDF. (B) The same histogram but with a 3-look PDF. 
The match is better. (C) The histogram of normalized amplitude ratio, HH/VV and 
the 3-look PDF. 



averaging pixels separated by a distance of two pixels. (A) HH-VV phase difference 
histogram and the 4-look PDF. The agreement is good. (B) The histogram of 
normalized |HH VV| and the 4-look PDF. (C) The histogram of normalized |HH|/|VV| 
and the 4-look PDF. 


REFERENCES 

[1] Goodman, J.W., 1985, Statistical Optics, John Wiley and Sons, New York, 1985. 

[2] Kong, J.A., et al., 1988, " Identification of Terrain Cover Using the Optimal 
Polarimetric Classifier," J. Electro. Waves Applic., 2(2) 171-194. 

[3] Lee, J.S., K.W. Koppel, S.A. Mango and A.R. Miller, 1993,"Intensity and Phase 
Statistics of Multi-look Polarimetric SAR Imagery," Proceedings of IGARSS’93. 

[4] Sarabandi, K., 1992, "Derivations of Phase Statistics from the Mueller Matrix", 
Radio Science, 27(5), 553-560. 

[5] Zebker, H. A., and Villasenor, J., 1 992, "Decorrelations in Interferometric Radar 
Echoes," IEEE Geoscience and Remote Sensing, vol.30, no. 5, 950-959. 


36 



